Quantum retrodiction in open systems 
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Quantum retrodiction involves finding the probabilities for various preparation events given a 
measurement event. This theory has been studied for some time but mainly as an interesting concept 
associated with time asymmetry in quantum mechanics. Recent interest in quantum communications 
and cryptography, however, has provided retrodiction with a potential practical application. For 
this purpose quantum retrodiction in open systems should be more relevant than in closed systems 
isolated from the environment. In this paper we study retrodiction in open systems and develop a 
general master equation for the backward time evolution of the measured state, which can be used 
for calculating preparation probabilities. We solve the master equation, by way of example, for the 
driven two- level atom coupled to the electromagnetic field. 
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I. INTRODUCTION 



The usual formulation of quantum mechanics is pre- 
dictive. It is well suited for calculating probabilities for 
particular outcomes of measurements from given prepa- 
ration information. The preparation information is nor- 
mally incorporated into a density operator representing 
the state of the prepared system, which evolves until the 
time of measurement. With sufficient knowledge we can 
assign a density operator in the form of a pure state 
projector to the system. This contains the maximum 
amount of information that nature allows us for predic- 
tion. The probability of the measurement outcome also 
depends on the operation of the measurement device. For 
a von Neumann measurement a particular measure- 
ment outcome is associated with a pure state projector; 
for a more general measurement the outcome is associ- 
ated with an element of a probability operator measure 
(POM) §. The latter are non-negative definite opera- 
tors that sum to the unit operator on the state space of 
the system being measured. The retrodictive formalism 
of quantum mechanics is used far less often than the pre- 
dictive formalism. Originally introduced by Aharonov et 
al. H in investigating the origin of the arrow of time, 
the retrodictive formalism involves assigning a state to 
the system based on knowledge of the measurement out- 
come. This state is assigned to the system just prior to 
the measurement and evolves backwards in time to the 
preparation event. The probability of a particular prepa- 
ration outcome can be calculated from a knowledge of the 
backwards evolved state and the operation of the prepa- 
ration device. While the results obtained arc fully con- 



sistent with predictive quantum mechanics plus inference 
based on Bayes' theorem [Q, the retrodictive approach is 
much more direct and provides a different insight into 
quantum mechanics. 
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Most applications of retrodiction involve closed sys- 
tems, where the time evolution is essentially unitary 
[||. An important emerging area in which retrodic- 
tive quantum mechanics will become increasingly impor- 
tant is in quantum communication [^|, including quan- 
tum cryptography |Q. Here a sender transmits a sys- 
tem to a recipient after preparing it in a particular state. 
The general communication problem is for the recipient 
to retrodict the prepared state from the outcome of a 
measurement and a knowledge of the operation of the 
preparation device. For closed systems, the evolution be- 
tween preparation and measurement is unitary and the 
intrinsic time symmetry significantly simplifies the prob- 
lem of calculating the retrodictive evolution. In practice, 
however, realistic systems will interact with an environ- 
ment into which information is irretrievably lost. This 
introduces extra uncertainty in addition to the intrin- 
sic quantum uncertainty associated with the preparation 
and measurement and removes the simple time symme- 
try associated with unitary evolution in a closed system. 
In open systems with a large environment the predictive 
evolution can often be described by a master equation. 
This describes information loss as the system propagates 
forward in time and becomes more entangled with an 
environment that is never measured. In this paper we 
investigate the derivation of a retrodictive master equa- 
tion H that describes irreversible loss of information as 
the system propagates backwards in time. As will be 
seen, this is not a simple time inverse of the predictive 
master equation. 
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II. PREPARATION AND MEASUREMENT 
WITHOUT TIME DELAY 

The probabilistic interpretation of quantum mechanics 
is usually expressed in terms of the theory of measure- 
ment. Retrodiction, however, relies on preparation being 
probabilistic so, in this section, we give a brief overview 
of the less familiar quantum theory of preparation and 
measurement. To dispense with time evolution at this 
stage, we consider an experiment in which one person, 
the preparer, prepares a quantum system and then an- 
other person, the measurer, immediately measures it. 
We assume that there are readout mechanisms on the 
preparation and measurement devices that indicate the 
preparation event and the measurement event that oc- 
cur. The preparer chooses which preparation event will 
occur and the measurer chooses whether or not to record 
the subsequent measurement event depending on what 
this measurement event is. If the measurement event j 
is recorded, then so too is the corresponding preparation 
event i. These events are recorded as a combined event 
(i, j). The experiment is repeated many times with the 
preparer choosing various states and a list of recorded 
combined events is made. The probability associated 
with a particular combined event on the list is defined in 
the usual way as the suitably normalized occurrence fre- 
quency. Both preparer and measurer have some control 
over the statistics, which contain elements of preselection 
and postselection. 

The connection with quantum mechanics can be in- 
corporated in the basic symmetric expression H for the 
probability that a combined event chosen at random on 
the list is (i, j): 



where 



and 



F Ar (i . ) = W 
Tr(Ar) 
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(1) 



(2) 



(3) 



Here Aj and Tj are non-negative definite operators as- 
sociated with the preparation and measurement events 
i and j respectively. They indicate the state in which 
the system is prepared or measured as well as any bias 
in preparation or in recording the events ||. The sums 
in (0) and (|^) are over respectively the measurement and 
preparation events that can be recorded. We say that the 
set of preparation device operators Aj describes mathe- 
matically the operation of the preparation device and the 
set of measurement device operators Tj describes the op- 
eration of the measurement device. Later, in equation 
(|l6|), we shall see that A, is P A (z)/S? lod where /5f led is 



the density operator for the prepared state and P A (i) 
is the a priori probability for this state to be prepared, 
which reflects any bias in the preparer's choice. Thus A 
is a density operator representing the best description we 
can give of the state in which the system is prepared if 
we know the operation of the preparation device but have 
no knowledge of the preparation or measurement events. 
Likewise T is a density operator representing the best 
description we can give of the state in which the system 
is measured if we know the operation of the measure- 
ment device but have no knowledge of the preparation or 
measurement events. 



A. Preparation and measurement probabilities 

The basic relation p) allows us to find expressions for 
various probabilities |9| in terms of preparation and mea- 
surement device operators Aj and Tj, We use super- 
scripts A and T to denote if the probabilities are based 
on knowledge of the operations of the preparation and 
measurement devices respectively. Where a probability 
is based on knowledge of a particular preparation event i 
or measurement event j, this is shown in the argument of 
the probability in the usual way by a vertical stroke pre- 
ceding i or j respectively. This stroke essentially means 
'if. A comma separating events means 'and'. Thus, for 
example, P Ar (j\i) is the probability that measurement 
event j is recorded if preparation event i is recorded, 
based on knowledge of the operations of the preparation 
and measurement devices. 

From the basic relation ([[]), with @ and (H|), we can 
calculate the following probabilities: 



P Ar (z) = ^P Ar ( l , J ) = 



Jr(A£) 
Tr(Af) 



Tr(Af) 



P- (j|z) 



P Al (hj) = Tr(A t r,-) 
P AT (i) Tr(Aif) 



p Ar m 



Tr(Af,-) 



(4) 



(5) 



(6) 



(7) 



Expression (m is the probability that, if an experiment 
chosen at random has a recorded combined event, this 
event includes preparation event i. Likewise (||) is the 
probability that the recorded combined event includes 
the measurement event j. Expression (||) is the probabil- 
ity that, if the recorded combined event includes event 
i, then it also includes event j. Expression (^) can be 
obtained by limiting the sample space to those events 
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containing i and is essentially Bayes' formula Like- 
wise (Q) is the probability that the preparation event in a 
recorded combined event was i if the measurement event 
is j. 

It is clear that we can multiply all the Aj by a con- 
stant value without changing the above probabilities and 
similarly for Tj . We make use of this flexibility to set 



and 



TrA = 1 



Trr = 1 



(8) 



(9) 



thus giving these operators the same normalization as 
density operators. 

Expression (Q) can be used for prediction. We can cal- 
culate the required probability if we know the preparation 
device operator (PDO) Aj associated with the prepara- 
tion event i and if we know the mathematical description 
of the operation of the measuring device, that is, every 
measurement device operator (MDO) Tj. Similarly we 
can use (0) for retrodiction from our knowledge of Tj and 
all the Aj of the preparation device. 

In the situation where no measurement is made on the 
system, there is only one possible measurement outcome 
and we can speak in terms of an effective measurement 
device whose operation is represented by a single 'no- 
information' MDO that is proportional to the unit oper- 
ator f§. Substituting this MDO for tj in (@) yields the 
retrodictive probability P which we write in this 

case as P A (i), where 



P A (i) = TrA, 



(10) 



Expression ( |To| ) , which depends only on the operation of 
the preparation device, is what we would normally call 
the a priori probability for the preparation event to be 
i, that is, the probability in the absence of knowledge of 
the preparation outcome, the measurement outcome and 
the operation of the measurement device. 



B. Biased devices 

Consider the special case where the operation of a 
measurement device is such that the probability P Ar (i) 
that a recorded preparation event is i in the absence of 
knowledge of the preparation or measurement outcomes 
is equal to the a priori probability P A (i) for i. Then, 
from 



and ©, 



Tr(Ajf) 
Tr(Af) 



TrA,- 



(11) 



If this is true for all possible preparation devices, that is 
for all Aj and A, then from (@) T must be proportional 
to the unit operator acting on the space of the system. 
We say that such measurement devices, which faithfully 



preserve the a priori probabilities of preparation events 
in the record, have an unbiased operation. Thus for an 
unbiased measuring device operation we can write 



Gl 



(12) 



where G is a positive constant and 1 is the unit operator. 
We then define 



n, = Tj/G 

which allows us to write (k5|) as 



P Ar (j\i)=Tr(pr% 



where 



'•pred 



= Aj/TrAj 



(13) 



(14) 



(15) 



This is just the usual, that is predictive, density oper- 
ator that we assign to the prepared state on the basis 
of knowledge of the preparation event. We note from 
(H) and (H) that flj are also non-negative definite and 
sum to the unit operator. Thus IT, form the elements 
of a probability operator measure (POM) || . Expression 
(fill) , which applies to unbiased measurement devices, is 
sometimes regarded as the fundamental postulate for the 
probabilistic interpretation of quantum mechanics Q . As 
we shall see later, even when the operation of the mea- 
surement device is biased, it is still convenient sometimes 
to use /5f Icd as defined by ([l5]) and in this case the predic- 
tive probability^) becomes (27) instead of (fuj). Equa- 
tions (|l^) and (|l0) give an expression for Aj in terms of 
the predictive density operator associated with prepara- 
tion event i and the a priori probability that this state 
is prepared: 



Aj = P A (i)p. 



prcd 



(16) 



This allows us to interpret A in as just the predictive 
density operator that we would assign to the prepared 
state if we knew the operation of the preparation de- 
vice but had no knowledge of the measurement or of the 
preparation outcome. 

Likewise, in the situation where we know only the sys- 
tem state space and are totally ignorant of the prepara- 
tion we are able to represent the operation of the prepa- 
ration device by a single 'no- information' PDO that is 
proportional to the unit operator. Substituting this PDO 
for Aj in M) yields 



P r (j) = T^ 



(17) 



We shall call expression (|l7|) , which depends only on the 
operation of the measurement device, the a priori prob- 
ability for the measurement event to be j, that is, the 
probability in the absence of knowledge of the measure- 
ment outcome and any preparation information. 

P Ar (j) is the probability that a recorded event is j 
if the operations of the preparation and measurement 
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devices are known but the actual preparation event is 
not known. From (|^) and (|lj) this will be equal to the 
a priori probability for j if 



Tr(ATj) 
Tr(Af) 



= Trr, 



(18) 



If this is true for all possible measurement devices, that is 
for all Tj and T , then A must be proportional to the unit 
operator acting on the space of the system. We say that 
such preparation devices, which faithfully preserve the a 
priori probabilities of measurement events in the record, 
have an unbiased operation. The retrodictive probability 
(0) reduces for an unbiased preparation device operation 
to 



P Ai = Tr(pf <'-,) 



(19) 



where Hj sum to the unit operator and are the prepara- 
tion POM elements, which are proportional to At, and 



/<'/" ' ./ ' ''' ./• 



(20) 



The symmetry between retrodictive and predictive ex- 
pressions is evident. It is worth remarking, however, that 
the usual experimental situation is asymmetric. This 
arises not from any intrinsic asymmetry in quantum me- 
chanics but from the asymmetry in the operations of nor- 
mal preparation and measurement devices. The opera- 
tion of the latter can usually be described by a POM and 
so ( |l4| ) is applicable; the former is usually not described 
by a POM so we must use the more general formula (0) 
and not (|l^) . This also results in an asymmetry between 
the forms of P AT (i) and P Ar (j Hn that P Ar (i) is just the 
a priori probability P A (i) in (lOh, which is independent 
of the operation of the measurement device, but P Ar (j) 
does depend on the operation of the preparation device. 
It is convenient, therefore, to express P Ar (j) in terms of 
P A (i) . The expression for this is 



p Ar (j) = Y, pAr (M pA ^- 



(21) 



Thus when the operation of the measuring device is un- 
biased, Bayes' theorem [|ll[, for example, can be written 
in terms of the a priori probability P A (i) as 



p Ar m 



P AL (j\*)P A (i) 
J2i P Ar {j\i)P A {i) ' 



(22) 



The concept of an unbiased operation of a preparation 
or measurement device is of particular relevance for prac- 
tical cases. For an unbiased operation of a measurement 
device, the occurrence frequency of the preparation event 
i in the list of combined events is proportional to the a 
priori probability for i. An unbiased measuring device 
allows the preparer to control the relative occurrence fre- 
quency of a preparation event. Likewise, an unbiased 
operation of a preparation device allows the occurrence 



frequency of the measurement event j to be proportional 
to the a priori probability for j, where this probability is 
based on knowledge of the operation of the measurement 
device but no knowledge of the actual measurement out- 
come or of the operation of the preparation device. An 
unbiased preparation device allows the measurer to con- 
trol the occurrence frequency of the measurement event. 
In practice the operations of many, but not all, measuring 
devices are unbiased but preparation devices generally 
have biased operations. This corresponds to preselection 
without postselection. For such cases, as outlined above, 
we can write the predictive probability for measurement 
outcome j given preparation event i in terms of a predic- 
tive density operator and an element of a probability op- 
erator measure (POM) for the measuring device as given 
by @. That is, for a measuring device with an unbiased 
operation, the preparation device operator A, converts to 
a density operator and the measurement device operator 
Tj converts to a POM element. We cannot in general, 
however, write the corresponding retrodictive probability 
as the inverse expression involving an element of a prob- 
ability operator measure for the preparation device and a 
retrodictive density operator. To do this we would need 
the operation of the preparation device to be unbiased, as 
occurs in some situations such as the Bennett-Brassard 
protocol for quantum key distribution (J, fz| |l2| . 



III. TIME EVOLUTION IN CLOSED SYSTEMS 

So far we have considered the preparation time t p to 
be the same as the measurement time t m and in (Q) we 
have used Aj, the preparation device operator at t p . To 
allow now for time evolution in a closed system between 
t p and t m we replace Aj and A in the probability (^) by 
Aj(t m ) and A(£ m ) where 



h(t m ) = U{t m ,t p )kiU ] '(t m ,t p ) 



(23) 



and a corresponding expression for A(t m ). Here U(t mi t p ) 
is the unitary time-shift operator. We might interpret 
this replacement as a modification of the original opera- 
tion of our preparation device to include an evolutionary 
period so that the new preparation device operators are 
Ai(t m ). The resulting probability expressions obtained 
from (Q) will be modified accordingly. Thus the predic- 
tive expression (Jq) becomes 



p Ar m 



tn 5 J 1 j 



Tt\U{t m ,t p )A4W{t m ,t p )T] 
and the retrodictive expression (Q) becomes 
TrlUit^QAiU^t^tp)^] 



p AL m 



Ti[U{t m ,t p )AW{t m ,t p )T 3 



(24) 



(25) 



From ([lfj) , noting that TrA^ is unchanged by a unitary 
transformation, we see (23) is equivalent to the usual 
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time development expression 

pf cd {t m ) = U(t m ,t p )pr d (t p )U\t m , t p ) (26) 

where f>Y ed (t p ) is the usual, that is predictive, density 
operator for the prepared state at t p . Using (|l5|) we can 
rewrite (^) as 



Tr[ P r d (t m )T] 



(27) 



A natural interpretation of (p7j) is in terms of the system 



being prepared in the state /5f red (£ p ) at t p and then evolv- 
ing forwards in time to become pf lcd (i m ) at the measure- 
ment time. 

At first sight, it may appear that the evolution has 
destroyed the symmetry between (^) and (^) which have 
now become (^4|) and (|2Jj). However by using the cyclic 
property of the trace we rewrite both (|24|) and ( p5|) in 
terms of Tj(t p ) where 



fj(t) = U\t m ,t)f 3 U{t m ,t), 



(28) 



with t p < t < t T 



formula (25) becomes 



Then, for example, the retrodictive 



p Ar m 



Tr[A,f 3 fa)] 
TrlAf,^)] ' 



(29) 



Expressions ( |28| ) and ( |29| ) can be interpreted in terms of 
a new operation of the measurement device which incor- 
porates a prior time delay as part of the measurement 
process. This operation is described by the set of mea- 
surement device operators Tj(t p ). 

We can use (§8|) and @ to write (j29j) as 



P Ar (z|j) 



Tr[A t/ 5frfa)] 
Tr[A,5f%)] 



(30) 



where 



-rctr 

Pi 



(t p ) = U^{t m ,t P )pf tT {t m )U{t m ,t p ). (31) 



Expressions (j|(]) and (|T]) can be interpreted as the sys- 
tem being assigned a retrodictive state p™ tr (t m ) at time 
t rn on the basis of the recorded measurement outcome 
and evolving backwards in time to the preparation time 
tp. 

Any of the above four interpretations of the process of 
preparation and measurement can be applied to both the 
predictive and retrodictive probabilities ( 24 ) and (^5j) . It 
is not necessary, for example, to use the predictive den- 
sity operator for prediction, nor the retrodictive density 
operator for retrodiction. Use of the forward-evolved pre- 
dictive density operator is convenient, however, when we 
know the outcome of a particular preparation and wish 
to retain flexibility in calculating probabilities for a range 
of possible measurement events with different measuring 



devices. Likewise, if we know the measurement outcome 
and wish to have flexibility in calculating probabilities 
for a range of possible preparation events with differ- 
ent preparation devices, it is convenient to calculate the 
backward-evolved retrodictive density operator. This is 
the purpose of this paper, except instead of finding the 
unitary evolution for a closed system, we wish to study 
the non-unitary evolution of a more realistic open system. 



IV. OPEN SYSTEMS 

In a preliminary paper |^] we outlined the derivation 
of an open system retrodictive master equation for the 
particular case where the operation of the measurement 
device can be described by a POM. In this section we 
examine in detail the more general case where the set of 
(renormalized) measurement device operators may not 
necessarily form a POM and also give a proof of the non- 
negativity of the retrodictive density operator at all times 
between measurement and preparation. Later we solve 
the master equation for a driven two-level atom. 

We are interested in a quantum system S interacting 
with an environment E. The initial state of the environ- 
ment at time t p is known but the state of the environment 
subsequent to this is never measured. Before examining 
the retrodictive problem we outline first the usual pre- 
dictive problem. In the predictive problem, in addition 
to the initial state of E, we also know the initial state of 
S and wish to predict the probabilities for the outcomes 
j of possible later measurements of S. For this we can 
use expressions ( p7| ) and (^6|) where /5f rcd (i p ) is now the 
predictive density operator for the combined system of 
E and which will be the product of the initial density 
operators for E and S: 



-prcd 



(tp) 



-prcd 

Pe 



-prcd 

[ Pi,s ■ 



(32) 



This initial state will evolve to a state which is in general 
entangled in accord with (|2^) where the unitary operator 
will now involve the Hamiltonian of the combined system. 
Each measurement device operator will be a product an 
operator Tj t s acting on S and an operator acting on the 
environment. As discussed earlier, the measurement de- 
vice operator for a non-measurement is proportional to 
the unit operator as this is effectively an unbiased mea- 
surement with only one possible outcome ||. Thus the 
predictive probability (p7|) becomes for this case 



p Ar m 



Tr E slpr d (tm)fs®iE} ' 



(33) 



where is the unit operator on the space of the the 
environment and T$ is the sum of Tj,s- The trace is over 
the environment and the system S states. By defining a 
reduced predictive density operator for the system S at 
the time t as Ifl3| 



Pff(t) = ^ E [U{t,tp)p^ a ® pr s a W{t,tp)] (34) 



^prcd 



pred/yt/ 
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we can write the predictive probability in terms of oper- 
ators acting on the state space of S. Substituting (26) 
and (53) with t — t m into (R3h gives 



p av (M = 



< pred (i, 



- j,S\ 



Trs[p!f(t m )T s } 



(35) 



which is of the same form as the closed system formula 
(B7T). Thus in order to find the predictive probability 
P Al (j\i) we need only calculate the evolution of the re- 
duced predictive density operator rather than the evo- 
lution of the combined system plus environment. If the 
environment E has a large number of degrees of freedom 
and is little changed by the coupling to S, then if the 
environment is not measured we can make use of the ap- 
proximation that, for any time t between preparation and 
measurement 13 



U(t,t P )p E red 



■ prcd U\t, t p ) 



"Pis 



-prcd 

Pe 



• prod (i). 



(36) 



The Markov approximation 1 13 then leads to a general 
Markovian master equation for p pr g d (t) that, in standard 
Lindblad form, is given by 



kPred 
Pi,S W 



ih-'iHsJf^t)] 

-^A q pYf{t)A\ 
q 



-Pis 



prcd (t)AlA q 



(37) 



where H$ is the Hamiltonian for the system S without the 
environment and A q is a system operator. This general 
form of equation incorporates the Markovian requirement 
and conserves both the non-negative definitcness and the 
trace of the reduced density operator |p^| . 

The expression for the retrodictive probability (p^) is 



p Ar m 



Tr ES [ki,i 



>pT%(t P )} 



Tr ss[As 



-prcd-A 

•Pe r ; 



(t P )] 



(38) 



where the backward-evolved measurement device opera- 
tor tj(t) is given by (p8|), As is the sum of A iy s and we 
have set the proportionality constant between the envi- 
ronmental preparation device operator and p^ cd to unity 
because this is the a priori probability that the environ- 
ment is prepared in this state. For our case, where the 
environment is not measured, the measurement device 
operator Tj will be proportional to f^s ®1 E - As Tj(t) 
evolves backwards in time towards t p , it changes from 
being factorizable at t m to become in general more en- 
tangled. 

We seek to express P Ar (i|j) in terms of a reduced retro- 
dictive density operator so that it has a form similar to 
©, just as P AT {j\i) in © is similar to ©. The sim- 
plest way to do this is first to define a reduced backward 
evolved measurement device operator and then to nor- 
malize this so that the trace is unity. From (|38|) we see 



that we need a reduced backward-evolved measurement 
device operator in the form 

t jt3 (t) cx Tr E [p^ ed U^(t m ,t)f j} s ® l E U(t m ,t)]. (39) 

We can now define a reduced retrodictive density opera- 
tor for the system S at the time t as 



r 3 -,s(*) 



(40) 



Then by substituting (|40|) and (|39|) into (|38|) we can write 
the retrodictive probability in terms of operators acting 
on the state space of S: 



p AL m 



Trs [Aj 



spfditp)} 



Tr s [A s p$(t P )] 



(41) 



Thus, if we know the reduced retrodictive density oper- 
ator and the operation of the preparation device, we can 
calculate the retrodictive preparation probabilities by a 
formula of the same form ( |30| ) as for a closed system. 
For this, we wish to find an appropriate master equation 
governing the evolution backwards in time. We should 
note that the reduced predictive and retrodictive den- 
sity operators in (|J) and (|40|) are not simply related 
to each other and thus the master equation for the re- 
duced retrodictive density operator is not immediately 
derivable from ( p7[ ) , for example simply by reversing the 
sign of the time. This, of course, is not due to some in- 
trinsic time asymmetry in quantum mechanics; it arises 
from the asymmetry in the environmental boundary con- 
ditions. For both the predictive and retrodictive cases we 
are assuming that the environment is prepared in some 
particular state but not measured. To obtain symmetry 
we would have to assume in the retrodictive case that the 
environment is measured to be in some state and that we 
have no information about its preparation. 

The most straightforward way of finding a retrodictive 
master equation for p™s r (i) is to find an equation for 

the backwards evolution of Tj t g(t) and then us e (Eol) to 
obtain the master equation. To do this, we use Q26|), the 
group property of the time-shift operator 



U(t m ,t p ) = U(t m ,t)U(t,t p ), 



(42) 



where t is a time between t m and t p , and the cyclic prop- 
erty of the trace to rewrite the numerator of the right 
side of @) as 

^Es[U(t, t p )p^ Cd ® /5^g d f)^(£, t p )U\t m , t)Tj.s 

®i E U{t m ,t)). (43) 

We leave the denominator of the right side of (|3^) as it is. 
Clearly P AT (j\i) does not vary as we change t between 
t p and t m and so the derivative of the numerator ( |4^ ) 
with respect to t must vanish. Using the weak-coupling 
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approximation ( pq ) and ( |39| ) we can write fl4q ) as pro- 
portional to 

Tr 8 \fif(t)t jiS (t)]. (44) 

The vanishing of the derivative of this with respect to t 
yields 



prcd 



Tr s [pP; s ca (i)r„ s W] = -TrslKs (45) 

Substituting (^) into (|5|) gives, after application of the 
cyclic property of the trace, 

Tr s [pPf(t)f jtS (t)} = Trs^^i-^- 1 ^,^^*)] 



AlA q t jjS (t)}}). 



(46) 



This is true for all /Sfs d (i), so the evolution equation for 
the reduced measurement device operator is 

<? 

-t , s {t)A\A q ~A\A q T, hS {t)]. (47) 

We can then find the retrodictive master equation for 
PTs(t) ^ substituting (|47]) into the time derivative of 
(|40|). This gives 

-Yy2A\?${t)A q -pfi{t)A\A q 
i 

-A\A q p^{t)\ 

- 2pfi{t)Tv s {pfi{t) Y&M) (48) 



as the master equation for pf^it)- 

It is not difficult to see that ([48]) conserves the trace of 
/5™g r (t), but ( |4^ ) does not preserve the trace of Tj t s(t)- 
The price of this preservation is that while ( |f7|) is a lin- 
ear differential equation, the master equation (Eq) is more 
complicated. While (|4|) can be solved directly for some 
simple cases, the solution is not always obvious. In gen- 
eral, however, (ji|) is reducible to a linear equation by 
substituting a variable-trace operator 

B(t) = r U '(t)exp ^2 J* Tr{p retr (0£[AJ, A q ]}dA , 

(49) 

which gives a linear equation for Bit). Then p retI (t) can 
be recovered from the solution as B(t)/TrB(t). The lin- 
ear equation obtained, however, is the same as (M7|) so 



B(t) is just proportional to f j,s{t)- It is thus more con- 
venient in general to calculate retrodictive probabilities 
by solving pTj ) for Tj t s(t) and then using 



p AL m 



Tr s [k s f jiS (t p )] ' 



(50) 



which is derivable from (|39D,(|3q) and (|28|), instead of 
solving @ for p)°^{t) and using @). 

To give a physical interpretation of Tj^sit), we first 
show that this is a non-negative operator for all t between 
t p and t m as follows. Consider an operator a(t), acting on 
the state space of S, defined to obey the predictive master 
equation ( |37[ ) for times t between (and including) T and 
t m where T is some time such that t p < T < t m . Let <r(T) 
= \u){u\ where \u) is some arbitrary pure state of S. Then 
a {t) will have the form of a predictive density operator 
and so cr{t m ) will be non- negative definite. From (|3S|), 
Fj,s(tm) is just Tj t St which is also non-negative definite. 
Thus 



Trs[a(t m )T hS {t m )} > 0. 



(51) 



It is not difficult to show directly from (|37|) with a(t) 
in place of Pi T g d (t) and from (47) that the time deriva- 
tive of Tr s[&{t)Tj t s{t)] is zero for t between T and t m . 
Consequently 



Tr s [a(T)f,- s (T)] > 0. 



(52) 



Writing the left side of (|5^) as (u\Tj y s(T)\u) shows that 
Tj t s(T) is non-negative definite for all times T between 
tp and t m . 

The non- negativity of fVs'(i) has two consequences. 
Firstly it follows from ( [40| ) that pj e £(t) is also non- 
negative. This, combined with the trace of /5jg r (i) be- 
ing conserved as unity, means that ( |48| ) is a legitimate 
master equation. The second consequence is that Tj^s(t) 
can be considered to be measurement device operators 
associated with measurement events j. Furthermore, if 
Tj t s(t) sum to be proportional to the unit operator on the 
state space of S it follows from ( |47|) that this sum is also 
conserved. Thus if the set of operators Tj : s{t m ), which 
are proportional to Tj t s ® 1b, describes the operation of 
an unbiased measuring device, so too will the set of op- 
erators Tj^it) and thus also Tj t s(tm)- Consequently the 
open system considered here allows a similar interpre- 
tation as does the closed system equation (^9|) with the 
measurement, or collapse of the state, taking place imme- 
diately after the preparation time t p . It should be noted 
that in the common interpretation, in which the opera- 
tion of the measurement device is described by Tj t s(tm), 
the initial environment state is considered to be part of 
the description of the operation of the preparation device. 
For the new interpretation, with the measurement taking 
place at t p in accord with (|50|), the initial environment 
state is considered to be part of the description of the 
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operation of the measurement device. This is analogous 
to a homodyne detection system in which the initial state 
of the local oscillator is considered to be part of the oper- 
ation of the measurement device rather than part of the 
operation of the preparation device along with the signal 
field state. The conservation of non-negativity allows us, 
in fact, to interpret (^7j) and ( p9| ) in terms of the measure- 
ment or collapse taking place at any time t between t p 
and t m with non-negative preparation and measurement 
device operators defined appropriately. The physical in- 
terpretation of the invariance of P Ar (j\i) under changes 
of t is that measurable probabilities are independent of 
when we choose the collapse time. This underlines the 
somewhat arbitrary nature of this concept. 

By comparing the predictive and retrodictive master 
equations (|37]) and (|4^), we see that when the interac- 
tion with the environment can be ignored, for example 
for short enough time intervals, the evolution is given 
by just the first term in each. The two equations for 
this unitary evolution are then the same. It should be 
remembered, however, that we have defined the time t 
so that it increases in both cases. For retrodiction, it is 
more natural to define a premeasurement time as 



so that 



T = U 



t. 



■ retr . , 

-p,, (*)• 



(53) 



(54) 



The unitary part of the retrodictive master equation in 
terms of r is then seen explicitly to be the time inverse 
of the corresponding part of the predictive master equa- 
tion. That the complete retrodictive master equation is 
not simply the time inverse of the predictive equation 
is illustrated by the following. If there is an extremely 
long time between preparation and measurement, both 
pP™ (t) and /5J e * r (T) can reach their steady-state values 

for which p[ g (t) and dp^ ct J (t) / dr are zero. A solution of 
the resulting steady-state retrodictive equation is clearly 
/5J e * r (oo) = ls/D where Is is the unit operator acting 
on the state space of S and D is the dimension of that 
space. This represents the no-information state. Essen- 
tially, at the measurement time the retrodictive state is 
the product of the measured state and the no- information 
environment state. As we go backwards in time from 
the measurement, the retrodictive system state becomes 
more entangled with the environment and we lose infor- 
mation about the system state. In the limit of very long 
times in the past, the system state becomes completely 
unretrodictable. There is no similar simple general so- 
lution of the predictive master equation, however, which 
would imply unpredictability for long times. Indeed for 
the case of an excited two-level atom undergoing spon- 
taneous emission into the environmental vacuum field, 
for example, it becomes very likely that the atom will 
be found in its ground state in the long term future, an 
outcome which is very predictable. 



V. COHERENTLY DRIVEN ATOM 
A. Retrodictive density operator 

As an example of a retrodictive master equation and 
its solution, we examine the case of a two-level atom 
driven on resonance by a strong laser field. The solu- 
tion of the predictive master equation for this system is 
well known fll3| , with the atom exhibiting damped Rabi 
oscillations at a frequency fl = (V 2 — 7 2 /4) 1 / 2 where V 
is proportional to the strength of the laser field and 7 is 
the spontaneous decay rate. We can write the Hamilto- 
nian describing the interaction between the atom and the 
laser field in the interaction picture in the semiclassical 
form 



(55) 



where a+ = \e)(g\ and tr_ = |<?)(e|. A comparison of (|37| ) 
with the predictive master equation for this system p3[ , 
that is, 



Pi,S (*) = 



iV 
— [(*+ 



-7[2*_^(i)S4 



a+a-Pls 



; prod (i) 



-prcd 

Pis 



(56) 



shows that q in ( p7|) has only one value, with A q — 
7 1 / 2 cr_. Thus, using B(t) in place of T^s{t), we can 
write the linear form (^) of the retrodictive master equa- 
tion in terms of the premeasurement time r as 

dB(r)/dr = l -^[(a + +a.),B(r)] 

-a + a.B(r)}. (57) 

We can solve ( p7| ) by converting it to a set of c- number 
equations. We write 



B(t) = u(t)&i + i>(t)<7 2 + w(t)<j 3 + x(r)lc 



(58) 



where Is is the unit operator acting on the space of the 
atom and 



<J\ = <J+ + <J- 

&2 = — i(d-+ — <7_) 
(T 3 = 2(7 + (T_ - 1 S . 

Then substitution of @ into @ yields 

du „ dv _ dw „ dx ~ 
T&1 + ~r^2 + —r^z + -r^s = 

dr dr dr dr 

— 7UCT1 — (jv — Vw)&2 — (Vv + 2~fw)&3 — 27wl^62) 



(59) 



(60) 



(61) 
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Multiplying (|6^) by o\ , a% and 03 respectively and taking 
the trace of both sides gives 



du 
rf7 



-7U 



dv 

— = —71' + Vw 
dr 



dw 

—— = — Vv — 271/7 
dr 



and taking the trace of both sides of (|6^) gives 
dx 

—— = —2jw. 
dr 

The solution of (|6^) is simply 

u(t) = u(0) exp(— jt). 



(63) 



(64) 



(65) 



(66) 



(67) 



The simultaneous equations ( p4| ) and[35|) are straight- 
forwardly solved by standard means |15|to yield 



v(t) = cxp(-37r/2){v(0)[cos(fiT) 
+ j(2Q)- 1 sm(QT)} 
+ u>(0)FfT 1 sin(fJT)} 



method involving solutions of the predictive master equa- 
tion and specific measurement POM elements. 

As our first example, suppose the atom is detected at 
time r = in the excited state |e)(e|. The measurement 
device operator corresponding to this detection event is 
proportional to |e)(e|, which is also the retrodictive den- 
sity operator at this time. Thus 5(0) oc |e)(e|. We note 
that it is not necessary to say whether or not the mea- 
surement device operator is an element of a POM. By 
writing |e)(e| as (1 + <r^)/2 we see that u(0) and v(0) are 
both zero and w(0) = x(0). With these values we find 
that 



u(t) = 



(72) 



v(t) = w(0) exp(-3 7 T/2)ra" 1 sin(fir) (73) 



W ( T ) = w(0) exp(-37r/2)[cos(fir) - 7(2fi)~ 1 sin(f2r)] 

* (74) 

and 



x(t) 



">(0) 
2 7 2 + V 2 



{V 2 - 2 7 exp(-37r/2) 



(68) 



h 2 + 2V 2 ) 
x[ 17 ^ in(»r)- 7 cos(Or)]}. (75) 



and 



w(t) 



exp(-37r/2){-w(0)ra" 1 sm(ftr) 
+ w(0)[cos(^t) - 7(2fi)~ 1 sin(fir)]} (69) 



which allows us to find x from (66) as 
2-y[Vv(0) - 710(0)] 



x(t) = x(0) + 



+ exp(-37r/2) 



2 7 2 + V 2 

■2 7 



2 7 2 + V 2 



x{ 



( 7 2 + 2V 2 )w{0) + 3jVv(0) 

2n 



sin(Or) 



-[Vv(0) -710(0)] cos(rJr)}. 



(70) 



Substituting (|67D, (|68|), Qfig) and (]70|) into fl5j) gives 
B(t). The retrodictive density operator is then just 



B{t) 
2x(t) ' 

which has unit trace as required. 



(71) 



B. Some detection events 

As a check of the retrodictive density operator (fn]) ob- 
tained from our general retrodictive master equation, we 
now calculate its matrix elements for some specific de- 
tection events. These can be compared with correspond- 
ing matrix elements in [ f[6[ calculated by a quite different 



The matrix elements of /5 ro * r (r) are easily found from ([H 
and (|5gh to be 



3,5 2x(t) 

w(t) + x{t) 
2x{t) 



(76) 
(77) 



with {g\p^{r)\e) = (e\p^(r)\g)* and {g\p\ et J {t)\ 9 ) 



1 - (e|pj e * r (r)|e). From @, @ and @ we obtain 



W:t r {r)\g) = exp(-3 7 r/2) sin(ttr) (78) 



2x{r)tt 



W^{r)\e) = 



w(0) 



■{V 2 



2x(t)(V 2 + 2j 2 , 
+ exp(-37r/2)[(y 2 + 4 7 2 ) cos(ftr 

7 



" ^(5^ 2 + 4 7 2 )sin(f7r)]}. 



(79) 



These are identical to the corresponding expressions de- 
rived in []l6| except that here we have x(t)/w(Q) in place 
of the normalization factor N. As we already have an 
expression for this in ( |75| ) there is no need to calculate it 
separately as we had to in fl(|. These results therefore 
not only confirm the validity of our new method but also 
display its advantages. 

The Rabi oscillations are clearly evident in both the 
diagonal and off-diagonal elements of the density matrix. 
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In the limit of long r, that is, in the distant past, the 
off-diagonal elements vanish and both diagonal elements 
tend to 1/2. This is the density matrix describing the no- 
information state discussed earlier. The actual state the 
atom was prepared in is thus essentially unretrodictable 
for these long times. This contrasts with the predictive 
density matrix, which also exhibits Rabi oscillations for 
short times but for long times in the future tends to an 
equilibrium state which is not the no-information state, 
instead it is determined by the relative value of V and 7 
& 

If the atom is detected in its ground state, we simply 
take u(0) and v(Q) as zero and w(0) = — x(0) and obtain a 
slight variation on the above elements of the retrodictive 
density matrix. 

If the atom is detected in the superposition state 
2~ 1 / 2 (|e) + \g)), the measurement device operator is pro- 
portional to (1 + <7i), so v(0) and w(0) are zero and u(0) 
= x(0). Then we obtain simply 



u(r) 
v(t) 
w{t) 
x(t) 



u(0) exp(— 71-) 





u(0). 



(80) 
(81) 
(82) 
(83) 



The retrodictive density operator is then easily found 
from (fnl) and @ to be 



p j re f(r) = [l + ( r 1 cxp(- 7 r)]/2. 



This has matrix elements 



(e\^(T)\g) = [exp(- 7 r)]/2 
<e|^ r (r)|e> = 1/2 



(84) 

(85) 
(86) 



which are precisely those calculated in |fLq| . Because of 
the phase of the detected state there are no Rabi oscil- 
lations, the density operator simply decays to the 
no-information state in the infinite past. To obtain Rabi 
oscillations in retrodiction from a detected state that is 
an equal superposition of \e) and \g), we could instead use 
a measuring device that detects the state 2~ 1 ' 2 (|e)+i|<?)), 
for which the measurement device operator is (1 + a 2) /2 
so v(0) is not zero. Some of the oscillatory terms in (68), 
( |69| ) and ((7^) are then retained. 



C. Preparation probabilities 

The elements of the retrodictive density matrix, di- 
agonal or off-diagonal, do not translate into preparation 
probabilities until we specify the operation of the prepa- 
ration device, that is, the values of A* £ for use with (41 ). 
The simplest calculations of preparation probabilities are 
those in the infinite-r limit. As we have seen, in this limit 
the retrodictive density operator becomes proportional to 
the unit operator. Substituting this into ( |4T| ) yields 



p Ar m 



Tr g A,,g 
Tr s A 5 ' 



which from (|§|) and (|1C|) is seen to be just the a priori 
probability for the preparation event i. Thus no further 
information has been gained from knowing the measure- 
ment event. 

We consider a preparation device which can prepare 
the two-level atom in a pure state \i) that is either |e), 
\g) or some superposition of the two. One might envisage 
a device that can apply pulses of resonant coherent radi- 
ation of various durations and phases to an atom in its 
ground state. The application of a 7r pulse would prepare 
the atom its the excited state and a 7r/2 pulse with ap- 
propriate phase would prepare the state 2~ 1 / 2 (|e) + \g)). 
A zero pulse would prepare the atom in the ground state. 
The operation of such a device is described by a set of 
preparation device operators, each of which represents 
a possible prepared state and incorporates the a priori 
probability for that state to be prepared. That is A^g is 
proportional to P A (i)\i)(i\. Clearly, even if the prepared 
states are restricted to |e) and \g) the preparation device 
operators will not be the elements of a POM unless these 
are produced with equal a priori probabilities. If the de- 
vice can only prepare the states 2 _1 / 2 (|e) + \g)) and \g), 
the preparation device operators are not proportional to 
the elements of a POM, whatever each a priori probabil- 
ity is. This is the situation for most preparation devices 
and can be contrasted with measurement devices, very 
many of whose operations can be described by a POM. 

To be specific, let us assume that the preparation de- 
vice is unbiased, preparing the atom in states |e) and \g) 
with equal probability. Then A e 5 = |e)(e|/2 and A S)1 g = 
\g){g\/2 and the preparation probabilities for states |e) 
and \g) become, from (p4[), simply the diagonal elements 
( e IPj°s r ( T )l e ) an d {g\p™f( T )\g)- On the other hand, for a 
biased device that prepares the atom in states |e) and \g) 
with a priori probabilities p and 1 — p respectively, A ei s 
and A 9t s are p|e)(e| and (1 — p)\g)(g\. Then the proba- 
bility that the atom was prepared in the excited state if 
the measurement is j is 



P AL (e\j) 



P(e\pT/(T)\e) 



p(e\p*f{T)\e) + {l-p)(g\p-f(r)\g) 



(87) 



In the limit of long r, this expression becomes equal to 
p, the a priori probability that the atom was prepared in 
the excited state, as discussed above. For short times r, 
the dependence of the probability (|8S| ) on r will exhibit 
an oscillatory behavior depending on the state detected. 

A simple but interesting case is for an unbiased prepa- 
ration device that prepares the states |+) = 2~ 1//2 (|e) + 
\g)) and |-) = 2^ 1 / 2 ( | e> — 1^> ) with equal probability. The 
preparation device operators are |+)(+|/2 and |— )(— 1/2. 
The preparation probabilities for +) and |— ) if the mea- 
surement event is j are, from (^Tj), just {+\p cc f and 
(— |/5J e * r (r)|— ) respectively. For the case discussed earlier 
in which the atom is detected in the superposition state 
2~ 1 / 2 (|e) + |oV) and so the retrodictive density operator 
is given by (B#, we find that these probabilities are just 
[1 + cxp(— jt)]/2 and [1 — exp(— jt)]/2 respectively. The 



11 



values of the probabilities change from unity and zero 
for small r to both being one half for large r and do not 
exhibit oscillations at any time. Thus the decay from 
the measured state to the no-information state is shown 
explicitly as the only change with time involved. 

The above examples illustrate the procedure involved 
and the preparation probabilities for preparation devices 
with other operations are easily calculated from (^l|). 

VI. CONCLUSION 

The emerging importance of quantum communication 
and quantum cryptography has made it worthwhile to 
re-investigate quantum retrodiction as a means of solv- 
ing the basic quantum communication problem. This in- 
volves calculating the probabilities that particular states 
were prepared when various states are detected. Where 
original investigations of quantum retrodiction were in 
terms of closed systems, quantum communication is more 
likely to involve open systems because of interaction with 
the environment. For closed systems the evolution is 
unitary, so the backward-time evolution equation for the 
state from measurement to preparation is just the simple 
inverse of the forward-time equation. For open systems 
the situation is not so simple. When a system is weakly 
coupled to a large environment, the forward-time evolu- 
tion is given by a master equation based on knowledge of 
the prepared state of the system and the initial state of 
the environment. In the retrodictive situation, the mea- 
sured state of the system is known, but not the final state 
of the environment because this is not measured. Instead, 
the initial state of the environment is known. Thus the 
backward-time evolution equation will not be the simple 
inverse of the forward-time master equation. In this pa- 
per we have derived the general, that is Lindblad, form of 
the backward-time, or retrodictive, master equation. We 
also prove that the retrodictive density operator remains 
non-negative definite at all times between measurement 
and preparation. As well as confirming the legitimacy 
of the master equation as an equation for the a retrod- 
ictive density operator, the non-negativity also allows us 
to interpret the evolution plus the measurement in terms 
of the operation of another measurement device. This 
shows that we can consider the measurement, or collapse 



of the state, as taking place immediately after the prepa- 
ration if we wish. Indeed it is possible to have a consis- 
tent interpretation of the preparation-measurement pro- 
cess with the collapse taking place at any time between 
preparation and measurement, showing the somewhat ar- 
bitrary nature of the concept. The vanishing of the time 
derivative of (pj), or in general of (|43|), effectively in- 
corporates the invariance of the physically measurable 
probability P Ar (j\i) with choice of collapse time. 

The retrodictive master equation is not in general lin- 
ear, underlining the fact that it is not the simple inverse 
of the predictive master equation. We have shown, how- 
ever, that it can always be linearized by means of a suit- 
able substitution, allowing analytical solutions for par- 
ticular systems. As an explicit example, we have consid- 
ered retrodiction for an atom driven by a coherent optical 
field and have solved the master equation for the general 
retrodictive density operator at all times between mea- 
surement and preparation. This allows us to select any 
particular measurement event we wish and easily find 
the probabilities of various possible preparation events 
based on this knowledge. Where the results overlap with 
previous calculations based on an entirely different, and 
less direct, method p6| there is perfect agreement, which 
serves as confirmation of the validity of our new ap- 
proach. 

As we have discussed previously m , preparation prob- 
abilities can always be found from predictive evolution 
plus inference by means of Bayes' theorem. Retrodictive 
quantum mechanics, however, offers a more direct and 
often simpler means of finding these probabilities as well 
as giving a different insight into the interpretation of the 
quantum state. In this paper we have extended the the- 
ory of quantum retrodiction for open systems by deriving 
of a general retrodictive master equation. 
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